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Abstract. We study the neutral minima of two-Higgs doublet models, showing that these 
^ . potentials can have at least two such minima with different depths. We analyse the phenomenol- 

QQ I ogy of these minima for the several types of two-Higgs doublet potentials, where CP is explicitly 

On ■ broken, spontaneously broken or preserved. We discover that it is possible to have a neutral 

. minimum in these potentials where the masses of the known particles have their standard values, 

I with another deeper minimum where those same particles acquire different masses. 
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1 Introduction 



X 



The standard model (SM) of electroweak interactions is a remarkably successful theory, but 
its scalar sector is yet untested. Numerous theories with a larger scalar content have been 
' proposed over the years in the attempt to increase the predictive power of the model and offer 

^ explanations for problems such as baryogenesis, CP violation, the hierarchy question and others. 

In models with more than one scalar, the possibility of the scalar potential having more than one 
minimum arises, and some of those minima can break the SU{3)c x C/(l)em gauge symmetry of 
the SM. Thus one is left with the possibility of imposing charge and/or colour breaking bounds 
on the parameters of the theory: excluding those combinations of parameters for which the 
potential's deeper minimum breaks charge and/or colour conservation. This appealing idea was 
first considered in the framework of supersymmetric theories [1] and applied extensively to such 
models [2j. Recently, it has also been applied to the Zee model [3] and to the two-Higgs doublet 
model (2HDM) [HIS]. The results obtained in this last reference were generalized in ref. [6j to 
models with an arbitrary number of Higgs doublets [^. 

The results of [4j may be summarised as follows: in the 2HDM, whenever a minimum that 
preserves charge and CP (we dub these "normal" minima) exists, that minimum is certainly 
deeper than any charge breaking (CB) stationary point. Further, the global minimum of the 
potential is a normal one, and the CB stationary point ends up being a saddle point. A similar 
conclusion holds for the spontaneous breaking of the CP symmetry: when a normal minimum 
exists, it is certainly deeper than any possible CP stationary point, and the global minimum is a 
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normal one. However, unlike the CB case, the question of whether the CP stationary point can 
be a minimum was left unanswered in ref. [1]. There is then the possibility that the potential 
may have a CP minimum, with a normal one lying below it. 

The importance of these results is to ensure the stability of the normal minimum against 
spontaneous charge or CP breaking: if one finds a normal minimum in the 2HDM, one may rest 
assured that there is no deeper CB or CP minimum for which the system may eventually evolve 
via tunneling. However, the work of ref. [1] did not answer the following question: how many 
normal minima does the 2HDM have? Are they all acceptable minima, with phenomenology 
according to the experimental data? And if there are several normal minima, which is the 
deepest one? Several of these questions were addressed for the 2HDM in recent papers [8], but 
in this paper we will focus on other aspects not treated in those works. A recent work [9j studied 
the minimum structure of the next-to-minimal Supersymmetric Standard Model. 

In the present work we will take a close look at the normal minimum structure of several 
types of 2HDM potentials. We will show that there is indeed the possibility that several normal 
minima coexist in the 2HDM, and if they are not degenerate, they will give rise to different 
phenomenologies. This paper is structured as follows: in section [2] we will introduce our formal- 
ism, write down the most general 2HDM potential and analyse the several types of theories one 
obtains by imposing symmetries on the model. In section [3] we will discuss the vacuum structure 
of the model and review the results of [3] about the differences in depths of the scalar potential 
at two of these possible vacua. In section |4] we will analyse normal minima in potentials where 
CP is broken - explicitly or spontaneously - and, in section [5l in potentials where CP is always 
left unbroken. We conclude, in section [6l with a general overview. 



2 The 2HDM potentials 

Let us consider two SU{2) doublets, with hypercharge Y = 1, given by 

<i>,= hl'.^') , <l>,= (^'l'.^A . (1) 

where all the (pi are real functions. Their numbering may seem odd, but it simplifies the writing 
of the scalar mass matrices. With these two fields one can build four SU (2) x U{\) real quadratic 
invariants, 

xi = l^ip =^1 + ^1 + ^1 + ^1^ 

X4 = /m($|^>2) = (fiipi - ip2^3 + - ^efr ■ (2) 

The most general 2HDM potential (for an overview, see for instance [?, iHj) is therefore a 
polynomial on the x's, with all possible linear and quadratic terms in these variables. That is, 

V = ai xi + 02 X2 + 03X3 + 04X4 + 611 xj + 622 + 633 x| + 644 x| + 

612 X1X2 + bi3 X1X3 + 614 X1X4 + 623 X2X3 + 624 X2X4 + 634 X3X4 , (3) 

where the coefficients Oj and bij are all real, the former having dimensions of mass squared and 
the latter being dimensionless. Under a CP transformation of the form <I>j — > xi, X2 and 
X3 remain the same but X4 switches signal. Thus, the terms of the potential which are linear 
in X4 break CP explicitly. The most general explicit CP-breaking potential has therefore 14 



real parameters. Through an appropriate choice of basis for {<I>i , $2} it is possible to reduce 
that number to 11 [12]. It is easy to see how to reduce the number of parameters to 12. The 
quadratic part of ([3]) may be written as 



aixi + 02X2 + 03X3 + 04X4 



$1 <I>^ 



ai 

03 + i 04 



03 — ta4 
2 

02 



^1 
$2 



(4) 



Through an unitary transformation on the fields, it is possible to diagonalize this 2x2 matrix. 
Written in terms of the new fields, the quadratic terms are now just two, whereas the quartic 
potential continues to have 10 parameters. In short, given the most general 2HDM potential, 
it is always possible to choose a field basis for which one has 03 = 04 = 0. To ensure that 
a minimum exists away from the origin, then at least one of the parameters ai, 02, has to 
be negative. Equation (fTSj) below shows us that one of the parameters {614 , 624 , ^34} may be 
expressed in terms of the two others, which leaves us with 11 independent real parameters. 

If we wish to study a potential for which CP is not explicitly broken, we must set to zero all 
terms linear in X4 in eq. ([3]): 04 = 614 = 624 = &34 = 0. Again, we may choose a basis in which 
03 = 0, and we are left with a 9-parameter potential. For this potential, one has the possibility 
of spontaneous breaking of CP. 

Finally, we may impose further symmetries on the model such that spontaneous CP breaking 
becomes impossible, and one is left with potentials containing seven independent real parameters. 
For instance, if one asks that the potential is invariant under a global Z2 symmetry ($1 —^1, 
$2 — > ^2), then the terms linear in X3 and X4 vanish and the potential becomes 



Va 



ai Xi + 02 X2 + 611 xl + &22 X2 + 633 xl + 644 xl + 612 X1X2 



(5) 



We call this the potential A. It is very important to have 633 7^ ^44 in this model, otherwise 
one has a massless axion0. Alternatively, with a U{1) symmetry (<I>i e*°$i, ^2 ^2), the 
model is also greatly simplified and one obtains the potential B, 



Vb = aixi + a2 X2 + 03 X3 + 611 x\ + 622 xl + ^33 {xl + X4) + 612 X1X2 , 



(6) 



where the symmetry imposes ^33 = 644. This potential is the analogous of the supersymmetric 
one. Notice that the 03 term softly breaks the U{\) symmetry. It was left there to prevent the 
appearance of a massless axion in the theory. Unlike the most general case, here one cannot 
"rotate away" the 03 term without introducing further parameters in the quartic terms of Ve- 
in order to ensure that a minimum exists away from the origin, a trivial calculation shows that 
we must have a\ > 4oi a2. As was shown in ref. [T3], the imposition of these symmetries makes 
it impossible for spontaneous CP breaking to occur in models A or B. Another property of 
these symmetries is that they are exactly those that are necessary to prevent flavour changing 
neutral currents (FCNC) in the theory, a phenomenon that plagues both the potential ([3]) and 
its restriction to the case of CP not being explicitly broken. 

Another possible symmetry to impose on the potential is that it remains invariant under an 
interchange between both flelds ($1 <-*• $2)- Then the potential becomes [14 



Vc 



ai {xi + X2) + 03x3 + 611 {xl + xl) + 633x1 + b4ixl + 

612X1X2 + 613 (Xl + X2)X3 . 



(7) 



Again, no basis changes can reduce the number of parameters in this model. 



^In fact, without this restriction, the potential has the (7(1) symmetry of model B without its soft breaking 
term. 



The potential Vc was first studied by Branco and Rebelo in [lAl, introducing an extra term 
of the form 02 (xi — X2) to softly break the symmetry imposed on the $ fields. The consequence 
is that the resulting potential may have spontaneous CP breaking, without FCNC. Similarly, 
adding to a term of the form 03X3, one breaks softly the Z2 symmetry that characterizes 
that potential. Once again, this gives rise to the possibility of spontaneous CP breaking. In 
fact, in ref. [12] it was shown that the potentials A and C are equivalent - it is possible to obtain 
the one from the other with a basis transformation. 

Interestingly, no soft breaking of the U{1) symmetry of the potential Vb gives rise to spon- 
taneous CP breaking. The potential written in eq. Q has already a soft breaking term - the 03 
term - and the only remaining term quadratic in the fields one could add to the potential would 
be 04 X4 - and that term explicitly breaks the CP symmetry. 

Let us summarise what we have discussed about the 2HDM potentials in terms of the number 
of independent real parameters that they have: 

• Potentials where CP is explicitly broken, which have 11 real parameters. FCNCs occur. 

• Potentials with 9 real parameters where explicit CP conservation has been imposed, but 
for which one might have spontaneous CP breaking. Again, FCNCs occur. 

• Potentials with 7 real parameters where, besides explicit CP conservation, one has also 
imposed a discrete symmetry {Z2 or permutation of <I>i and $2)- No spontaneous symmetry 
breaking occurs, FCNCs are excluded. With an additional term that breaks softly the 
discrete symmetry, spontaneous CP might arise. 

• Potentials with 6 real parameters, with explicit CP conservation and invariant under a 
global U{1) symmetry. Again there is no FCNC but these models have a massless axion. 
If one adds to the potential a term that softly breaks the global U{1), that axion acquires 
a mass but there is no possibility of spontaneous CP breaking. 

There is an interesting possible classification of these potentials, which we will present in ap- 
pendix [HI 



The vacuum structure of the 2IIDM is very simple |10j . with three different types of vev con- 
figuration. For a charge and CP conserving vacuum, the fields $1 and $2 have vevs of the 
form 



where vi and V2 are real numbers. For a vacuum that breaks charge conservation, we have 



Again, all vevs in this equation are real. Finally, a CP breaking vacuum occurs when the fields 
acquire vevs of the form 



3 The stationary points of the 2HDM 




(8) 




(9) 




(10) 
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Vacua with a and 6 simultaneously non-zero are not considered because the minimisation con- 
ditions of the potential forbid them H. 

The structure of the potential of eq. ([3]) is such that it may have several stationary points, 
and they may be of different natures. In ref. [1] we obtained a remarkable result that relates the 
difference in the depths of the potential in each of the three types of possible stationary points. 
Namely, the difference between the value of the potential at a CB stationary point, Vcb, and 
the value of potential at a normal stationary point, Vn, is given by 

Vcb - Vn = [iv[v2 - v'^vif + vl] , (11) 

where v"^ = v\ + f| and -/Vf^± is the value of the squared mass of the charged Higgs scalar, 
evaluated at the normal stationary point. What this equation tells us is that when the normal 
stationary point is a minimum H then, all of the eigenvalues of the squared scalar mass matrices 
being positive, we will have > and the deeper stationary point will be the normal 

minimum. Furthermore, in ref. [1] we were able to prove that the CB stationary point is unique, 
and that eq. (jlip . when > 0, also implies that it is, necessarily, a saddle point. The 

stability of the normal minimum against tunneling to a deeper charge breaking stationary point 
is thus ensured in the 2HDM. 

A similar result holds for the comparison between the CP and normal stationary points. In 
the case where it makes sense to define a CP symmetry - that is, in potentials where it is not 
explicitly broken - the difference between the value of the potential in the stationary points with 
vevs given by eqs. ([8]) and ([TO]) is given by [1| 

VcP - Vn = ^ [{v'lv2 - v'iv.f + 5^vl] . (12) 

M\ is now the value of the squared pseudoscalar mass at the normal stationary point. Again, if 
this squared mass is positive, for instance if the normal stationary point is a minimum, then the 
normal minimum is indeed deeper than the CP breaking one. Again, the CP stationary point 
being unique, the stability of the normal minimum against tunneling is guaranteed. However, 
unlike the CB case, we were not able to prove, in ref. that the CP stationary point is 
necessarily a saddle point when the normal stationary point is a minimum. Thus, even though 
we have proved that the global minimum is a normal one, we cannot discard the possibility that 
above it a CP breaking minimum exists. 

The final case of interest occurs in the potential with explicit CP breaking. In that situation, 
both the vevs of eq. ^ and those of eq. (jlOp have the same physical meaning: both stationary 
points break the same symmetry and the existence of a relative phase between the vevs of each 
field does not distinguish them. For this potential, we refer to the vevs of ([8]) as the 'Wi 
stationary point" , and those of (jlOp as the 'W2 stationary point" . It may be shown [12] that it 
is possible to pass from a vev structure of the type N2 to vevs of the A^i type by a specific basis 
transformation - the complex phase of eq. (jlOp is absorbed in the parameters of the potential, 
and the final vevs of both $1 and $2 are real. However, this does not mean that the A^i and 
stationary points are the same - simply that there is a field basis for which may be written 
with real vevs. In that new basis, however, the A^i vevs would have acquired a relative complex 
phase. In ref. [4J we found an interesting relationship between the difference in the depths of 

^Except for a very special case in the explicit CP breaking potential. Even in that case, though, via a basis 
change, that vacuum may be reduced to one with a 7^ and 5 = Q. 

^Even more generally, when it is a stationary point for which > 0. 



the potential at the A^i, ^^^^2 stationary points. Namely, 



— Vni 



Ml, 



Ml, 



N2 



[( 



In this equation we have {v'^)ni 



+ V2 and (u^ 



IN2 



vf + vf + 52 



and (M^ 



(13) 



are 



the squared charged scalar masses at each of the A'^i, stationary points. A similar expression 
was found in ref. [§]. Equation (jlSh tells us the deepest stationary point will be the one with 
the largest ratio between the square of the charged Higgs mass and the v'^. 

Let us now consider the stationarity equations that give rise to the different stationary points 
that we have been discussing. We have mentioned that the CB and CP stationary points are 
unique since they are given by linear equations on the vevs. However, this is not true for the 
normal stationary point. Let us begin with the most general 2HDM potential written in a basis 
where 03 = 04 = 0. We make ip^ = vi and (^e = V2 in eq. and minimise by solving 
the equations dV/dipi = 0, i = 1 . . . 8. Most of these equations are trivially satisfied. The 
non-trivial ones are 



dV 

dV 

dipQ 



2aivi + ibuvf + 2(612 + b33)viv'l + 36i3t'fv2 + 623 = 
2 02 V2 + 4 622 + 2 (612 + 633) vj V2 + 613 vf + 3 623 V1V2 = 



(14) 



and 



dV 

dip7 
dV 

dm 



vi (614 vf + 624 V2 + 634 Vi V2) 
- V2 (614 vl + 624 vl + 634 Vl V2) 



(15) 



Notice that one cannot have solutions of the form {vi = , f 2 7^ 0} or {vi 7^ , f 2 = 0}, 
unless some parameters of the potential are set to zero (623, 624 and 613, 614 respectively). 
Since there is no symmetry forcing those parameters to be zero, they have to be present in the 
potential. We now define the usual polar coordinates vi = v cos (3 and V2 = v sin (3. A trivial 
solution of these equations is clearly w = 0. Excluding that case, the stationarity conditions (fHj) 
become 

2 1 2 oi 



cos2/3 623 tan^/? + 2(612 + ^33) tan^/? + 3613 tan/3 + 4 611 



(16) 



and 



- a2 623 tan^ (5 + [4 oi 622 - 2 02 (612 + 633)] tan^ [3 

+ 3 (ai 623 - 02 613) tan^ (3 + [2 oi (612 + 633) - 4 02 611] tan /3 + oi 613 = (17) 



and both equations p5|) reduce to 

624 tan^ (3 + 634 tan (3 + 614 







(18) 



Eq. (jl6p tells us that, other than its sign, the value of v is determined unequivocally by tan j3. 
Eq. (fT7|) is a quartic equation on tan /3, having at most four possible real solutions. These two 
equations describe therefore eight possible solutions {vi , ^2}) due to the ambiguity on the sign 
of V. The 2HDM potential ([3]) is however invariant under the transformation <I>i — > — <I>i and 
$2 — ^ —"^2) so that these eight solutions correspond to only four different physical scenarios. 
Adding the trivial solution f 1 = 'L'2 = 0, we have a total of nine solutions. However, we must 
contend with eq. (jl8p as well, which is a quadratic equation on tan (3. Then, there are at most 



two different values of tan /3 wliich satisfy all equations. This means that we have a maximum 
of five stationary points. 

For potentials with explicit CP conservation, equations (jlSp are trivially satisfied, since 
6i4 = 624 = ^34 = 0. Therefore, equation (llSp does not exist and the potential could have a 
total of nine stationary points. 

At this point we ask: can we have more than one normal minimum, with different depths? 
The answer is yes, and to see this we make use of Morse's inequalities [15j: for a given real 
function of two variables, let niQ, mi and m2 be the number of its minima, saddle points and 
maxima, respectively. For a polynomial function in vi and V2, bounded from below, such as the 
one we are dealing with, Morse's inequalities state that: 

• mi > m-o — 1; 

• mo — mi + m2 = 1. 

We know that the 2HDM potential has mo + mi + m2 = 2n + 1 stationary solutions, n = 
0, . . . , 4: at most 2n real roots of eqs. (fT6|) . (fT7|) and (fT8|) plus the trivial solution vi = V2 = 0. 
Hence we find that mo + m2 = n + 1. Let us analyse the several possibilities for the number 
of minima mo, depending on the number of solutions n. Simply counting all the different 
combinations of extrema leads us to: 

• n = 0: we have necessarily mo = 1, the minimum is unique but is located at the origin, 
vi = V2 = 0, which means that there is no SU{2)y\/ x U{1)y symmetry breaking. This 
case is excluded on physical grounds. 

• n = 1: we find two possibilities: 

— mo = 1, which is the previous case. 

— mo = 2, which means two degenerate minima away from the origin, related to one 
another by a change of sign of the vevs. This situation corresponds to an acceptable 
symmetry breaking and it means that there are no normal minima with different 
depths. This would be the "standard" situation. 

• n = 2: there are three possibilities: 

— mo = 1 or 2 are like the previous cases. 

— mo = 3, one uninteresting minimum at the origin, and two degenerate ones away 
from it. This situation would also be the "standard" one, as there would be no normal 
minima with different depths. 

• n = 3: there are now three qualitatively different cases: 

— mo = 1,2,3, like above. 

— mo = 4, this case corresponds to two pairs of degenerate minima away from the 
origin. Nothing forces these two pairs of minima to have the same depth. We might 
therefore have one normal minimum deeper than another. 

• n = 4: we have: 

— mo = 1,2,3,4, like above. 



— mo = 5 is similar to the mo = 4 case examined above, with an extra minimum 
present at the origin. 



This trivial analysis shows us that, if there are more than two solutions for tan /3, then the 2HDM 
may have more than one normal minimum away from the origin at different depths. However, 
no more than two such minima can exist. 

Why is this interesting? We already know that in the 2HDM, when a normal minimum 
exists, then the global minimum of the theory is normal [4J. However, another interesting 
possibility might arise: the minimisation equations may give us a normal minimum A^i with vevs 
{vi , V2} and another normal minimum N2 with vevs {vi , ^2}- Which is the deepest? Suppose 
that we find A^i for which the vevs are such that the SM phenomenology is satisfied, namely, 
V1+V2 = (246 GeV)^ and the W and Z boson masses are according to their experimental 
values; but N2 is deeper, for which the sum of the squared vevs vf + has a completely 
different value, contrary to experimental data. This situation is clearly undesirable since it can 
lead to tunneling between an acceptable minimum and one with gauge boson and fermion masses 
different from their measured values. A trivial calculation (see appendix [A|) shows us that the 
difference in depth of the potential at A'^i and is given by 



a result very similar to eq. (jl3p . Again, it is the squared mass of the charged scalars divided 
by the respective v'^ which "controls" the difference in depths of the potential at the stationary 
points. As we will shortly show, it is possible, depending on the value of the parameters of 
the potential, to have such "coexisting" minima in situations of interest for particle physics 
phenomenology. 

4 Neutral minima in potentials with CP breaking 

As mentioned before, one of the most interesting features of the 2HDM is that it allows for 
breaking of the CP symmetry. We will now analyse the possibility of coexistence of neutral 
minima in this potential, for the two cases of CP breaking: explicit and spontaneous. 

4.1 Potential with explicit CP breaking 

As we saw earlier, for the 2HDM potential with explicit CP breaking the value of tan/3 is 
determined by a set of two equations, (fT7|) and ([18|), one quartic and another quadratic. This 
means that we may have, at most, two different values of tan /?. We are therefore in the n = 2 
case discussed above. As was shown, this means that the normal minimum exists and is unique. 

This conclusion, however, was derived for stationary points with real vevs, of the form ([8]). 
We know that in this potential CP is not defined and as such stationary points of the form of 
eq. (jlOp . in which the vevs have a relative complex phase, have the same physical relevance as 
stationary points with real vevs. Equation (jl3p gives us the difference in depth of the 2HDM 
potential with explicit CP breaking at two neutral stationary points, A'^i (with real vevs) and 
N2 (with vevs with a relative complex phase). This result was first obtained in ref. [4], but one 
question was left open in that work: is it possible to find values of the potential such that one 
can find two minima that verify eq. (jl3p ? We will perform a numerical study to prove that such 
a situation is indeed possible. 




(19) 



O 



The stationarity conditions for a A'^i stationary point with vevs such as those of eq. ([8]) 
were shown in eq. (jl4p . Similarly, for a A''2 stationary point, with vevs of the form (jlOp . the 
stationarity conditions are 



^33)'fl ^^2 



bi3 {3vf 



2ai< + Abuv'livf + 6^) + 2 (612 + 

- 2 614 v'l V2 6 - 634 vf 6 

2 02 4' + U22vf + 2 6i2^;^'K^ + + '^h^vfv'i + bi^v'i{vf + 5' 

+ 3623^1 '^2^ — &14 (^1 ^ + 5^)5 — 36241^2^'^ — 2 634 ^"^2 5 + 2 644t'2(! 

2ai(5 + 4 611 5 + + 2 6i2 5^;f + 2 6i3<z;^'(5 - bu{vf + 3(5^) 

- 024 ^^2 - O34 V1V2 +2 644 ^2 



I //3 

023 ^2 



. 

(20) 



Because we are looking for simultaneous A'^i and A'^2 minima, these equations have to be solved 
together with the stationarity conditions and ([15]). There are two ways to do this: (a), 
generate a set of {oj , hjk\ parameters and use both Ni and A'^2 stationarity conditions to deter- 
mine the 5 vevs; or (b) generate all but 6 of the parameters {oj , hji^} and also the 5 vevs and 
use the stationarity conditions to determine the remaining 6 potential parameters. We chose 
option (b) for two reasons. Firstly, we want one of the minima to describe the "real" world. As 
such, we know what the sum of the squared vevs should be ((246 GeV)^) at that minimum. It 
is easier to use this information by inputing the vevs to begin with. Secondly, the stationarity 
conditions (jl4p . (j20p are cubic in the vevs, and therefore very difficult to solve analytically. 
Numerical methods need to be used, which are difficult to control and time-consuming. Deter- 
mining the potential parameters requires nothing more elaborate than solving a set of linear 
equations, once the vevs have been specified. Our method to solve the stationarity equations, 
then, consists in: (1) generating values for the vevs, {vi , ^2} and {v'l , V2 , 5}, and all but 4 
of the bij] (2) using eqs. (fill) to determine the value of the parameters {ai , 02}; (3) choose 
the parameters {614 , 624 , 634 , 644} so that they satisfy the remaining stationarity equations, 
namely: 



2 v'l v'^ 5 
6ivf + 6^) 





3v"l5 

//3 
V"2 

^2 



//2 r 

2v'iv'i5 
v'{v"l 

Vi V2 





-2v"2 5^ 
-2 6v"l 



14 




'Si 


24 




S2 


34 




S3 


44 




S4 



(21) 



with 



^3 
^4 



Si = 2aiv" + 4biiv"{vf + 6^) + 2(612 + 633) t-i^ + 613 (3 + S)v'^ 



,.'/2 



+ b23vf 



S2 = 2a2v'i + U22vf + 2bi2v'i{vf + 5^) + 2 633<^< + 613 < (< + 5^ 



+ 3 623^^l^^f 

2ai8 + Abii5{vf + 5^) + 2bi2 5vf + 2 613 7;^' w^' 5 
. 



(22) 



We chose random values for the potential's parameters, such that the bij couplings were of the 
same order - we considered b parameters in the range 10~^ < bij < 10. To be certain that the 
solutions of the stationarity conditions correspond to minima, we calculated the eigenvalues of 
the squared scalar mass matrices and verified that, except for the three zeros corresponding to 
the Goldstone bosons, the remaining ones are positive. The mass matrix expressions for any 
2HDM potential may be found in ref. [4j. 



Let us see, for instance, if it is possible that the A^i and A''2 minima have the same value 
for the squared vevs, but different scalar masses. In other words, can the 2HDM accommodate 
two minima which predict the same W, Z and fermion masses, but different scalar spectra? 
The answer is, yes. To see this, we solve eqs. ([22]) by inputing values for the vevs such that 
vf + V2 = v'l^ + + (5^ = (246 GeV)^. We also input the several bij parameters and scan 
the parameter space, accepting only sets of parameter values for which both A^i and N2 are 
minima. To make these minima of some physical interest, we also demanded that all the scalar 
masses be larger than 100 GeV, but inferior to 1 TeV. The results we found are illustrated in 
fig. ([I]), where we plot the difference of potential depths at A^i and N2 against the difference 
in the charged scalar masses at both minima. One immediately observes that the A^i and N2 
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Figure 1: Difference in depth of the potential at two minima, A^i and for the 2HDM model 
with explicit CP breaking, in terms of the difference of charged scalar masses. 

minima are equally likely to be the deepest minima. We also remark that it is not difficult to 
find combinations of {a , 6} parameters and vevs for which these minima coexist: for 10000 sets 
of parameters for which A'^i is a minimum obeying the criteria described above, about 10% have 
a N2 minimum "alongside" . Hence, for the same set of parameters of the potential we can have 
two different minima, both of them predicting the same values for the gauge boson masses and 
fermions, but with different spectra of scalar particles. Tunneling from A'^i to would therefore 
only change the values of Higgs' masses. 
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4.2 Potentials with spontaneous CP breaking 



For the 2HDM where CP is not exphcitly broken, ah of the couphngs 6^4 are set to zero, except 
for 644, so that the stationarity condition of eq. (jlSp is trivially satisfied. Then the number of 
possible normal minima goes up - we have a total of nine possible stationary points from which, 
as explained above, we may have a maximum of two non-degenerate minima away from the 
origin. For these potentials, though, the stationary points of the type N2 now correspond to a 
spontaneous breaking of the CP symmetry. 

For illustrative purposes, let us see if the following situation might occur: a CP breaking 
minimum with v = 246 GeV above a normal one @. At this normal minimum, the value of 
V is not fixed a priori. To achieve this end we once again solve the stationarity conditions of 
the 2HDM potential by generating random values for most of the parameters and vevs and thus 
obtaining linear equations on the remaining unknowns. Solving the linear set of equations thereof 
resulting produces a complete set of parameters for the potential, for which we then proceed 
to investigate whether both stationary points - the normal and the CP one - are minima, by 
analysing the eigenvalues of the corresponding squared scalar mass matrices. 




Figure 2: Masses of the neutral Higgs bosons at a CP breaking minimum with v = 246 GeV 
versus the expected W mass at the deeper normal minimum. For each combination of potential 
parameters, the lightest Higgs mass is represented by crosses, the second lightest mass by stars, 
the heaviest one by circles. 

The results we found are shown in fig. ([2]), where we plot the masses of the neutral Higgs 

^Since, we repeat, according to the results of ref. ^ the reverse is not possible. 
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bosons at the CP minimum, versus the W boson mass expected at the normal minimum. For 
visuahsation purposes, the vertical axis was limited to values above 10~^ GeV, which removed 
from view several points corresponding to very light Higgs masses. The normal minimum is 
found to have depth similar to that of the CP one, but the resulting W masses are in general 
very different. An interesting feature of the CP violating minimum is also shown in this plot: the 
neutral Higgs masses are extremely low, much lower than the most recent experimental bounds, 
with the exception of the heaviest neutral mass, represented by circles in fig. ([2]), which may 
reach hundreds of GeV. The second heaviest neutral Higgs (stars in fig. ^) is also very light. 
For the lightest neutral Higgs, the maximum value we obtained for its mass was of the order of 
0.1 GeV. One must also emphasize that these sets of parameters where CP minima coexist with 
normal ones are extremely difficult to find. Out of a total of about 10 million normal minima, 
we only find CP ones with a sum of squared vevs equal to 246 GeV for some 160 of those. Our 
conclusion, then, is that though coexisting normal and CP minima are numerically possible, they 
are extremely unlikely, at least in the range of parameters we chose. When they are found, the 
CP minimum seems to be characterized by having a very low mass for the lightest Higgs boson. 
This would seem to exclude the strange possibility that we are "living" in a CP minimum with 
standard W masses, with a normal minimum lying underneath it. Notice, however, that this 
does not exclude, at all, the possibility that the potential has parameters such that the global 
minimum of the theory is CP violating. The conclusion of [1] only imply that in this case no 
normal minimum exists. 

We have also considered the situation where the normal minimum has v = 246 GeV, with 
a CP minimum above it, with vevs determined by the stationarity conditions. In cosmological 
terms this would raise the possibility of a universe "resting" for some time in such a CP breaking 
minimum immediately after the Big Bang, before tunneling the present normal minimum. Again 
we found that this situation only occurs for a very small subset of all the parameter space that 
was scanned. The results obtained are very similar to those shown in fig. ([2]). The W masses 
found at the CP minimum and the neutral Higgs masses at the normal minimum tend to be 
smaller than 100 GeV. Nevertheless, we have found several cases where the Higgs masses are 
compatible with present experimental bounds. 

5 Normal minima in potentials where CP is conserved 

In the previous sections we looked at the simultaneous existence of two neutral minima, one with 
real vevs, the other with complex ones. As we showed in section [Sj however, it is possible to have 
several normal minima with real vevs {vi , V2}, not all of them having the same depths. This is 
already a possibility for the potentials with soft CP breaking, analysed in the previous section, 
but the case we treated there was the most interesting situation of a CP-breaking minimum 
coexisting with a normal one. However, for potentials where CP breaking is impossible the 
existence of several normal minima for the same potential is of interest. As was explained in 
section [2l there are two different potentials where spontaneous CP breaking is not possible, 
depending on the symmetries one has imposed on the 2HDM model. As we will now show, the 
normal vacua structure is very different for these two potentials. 

5.1 Normal vacua in the potential Va 

The stationarity conditions for the potential Va_ may be obtained from the most general ones 
written in eqs. (|14p and (llSp if one makes 613 = ^23 = 0. Since for this potential 614 = 624 = 
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634 = 0, eq. (fT5|) is trivially satisfied and because 613 = 623 = 0, eqs. (fHll become 



dV 

— = 2aivi + Abnv'f + 2{bi2 + b33)viv^ =0 

= 202^2 + 4622^1 + 2(612 + 633)^1^^2 = . (23) 

We can identify the trivial solution, vi = V2 = 0, and three non-trivial ones: 
• Solution I (if 02 < ol): 

.1 = 0, vl = . (24) 

^ 022 

At this stationary point the value of the potential is 

-"—it- <^^> 



Solution II (if oi < 0): 



2611 

At this stationary point the value of the potential is 



= - ^ ' ^2 = . (26) 



Solution III: 



^2 _ 2b22ai - (&12 + 633) «2 ^2 _ 2611 a2 - (^12 + 633) Qi ^28) 
^ (&12 + 633)^ - 4611622 ' (612 + 633)2 - 4611622 

At this stationary point the value of the potential is 

y^^^ _ 62201 + 611 a| - (612 + 633)0102 ^^g^ 
4 611 622 - (612 + 633)2 

The solutions I and II correspond to models very similar to the SM although, for certain realiza- 
tions of the 2IIDM, these solutions may correspond to cases where either the up or down quarks 
are massless. The solution III is the most interesting one, and we will analyse it in detail. To do 
so, let us first consider the relations one can establish between the {o , 6} parameters and some 
of the physical parameters of the model at a stationary point III: the neutral CP-even scalar 
masses Mh and Mh, the angle /? with tan/3 = V2/V1 and the rotation angle a that diagonalises 
the matrix of the squared masses of the CP-even scalars. These relations are found in ref. [16] . 
and are given by 

oi = — [cos a cos(/3 — a) Mfj — sin a sin(/3 — a) M|] 

4r COS l) 

02= — : TsiuQ cos(d — a) M% + cos a sin(d — a) M?] 

4 Sm/3 \r- ; ti ; m 

^11 = -r^ 0-7T cos aMrr + sin a Mi, 

4 f ^ cos^ p ' 

b22= , o ^ [WaMl + cos^aMl) . (30) 
4 f ^ sm p 



^Requiring that the potential be bounded from below implies that 611 > and &22 > 0. 
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There are more relations, but these four will be all that we will need for our purposes. Let us 
assume that we have a set of parameters of the potential for which the solution III defined above 
exists and is a minimum, and for which the solutions I and II exist as well. The question we now 
ask ourselves is, under these circumstances, is it possible that either I or II are also minima, and 
deeper than III? Let us start by comparing the value of the potential at the stationary point 
I and the minimum III. Making use of the fact that at any stationary point the value of the 
potential is given by (oi vf + 02 V2)/2 and using eqs. (pOl) . we may rewrite Vj and Vm as 



Vi= - — 



v'^ [sin a cos(/? — a) + cos a sin(/3 — a) Mj 



212 
hi 



4 sin^ aM^ + cos^ a M| 

Vm = - — [cos2(/3 -a)Mfj + sin^ip - a) Ml] . (31) 

Notice that we are rewriting the values of the potential in terms of the physical parameters 
(Higgs masses, etc) of the minimum III. Can we have Vj < Vjjj? According to the expressions 
above that occurs if 

Tsina cos(/? — a) + cos a sin(/? — a) M?l ^ „^ ^ „ ^, , r, 

^ ^. 2 J2 ^ > cos2(/3-a)Mi + sin2(/3-a)M2 ^ 

sm a Mjj + cos^ a 



[sin Q cos(/3 — a) + cos a sin(/3 — a) M|] > 

> [cos^{P -a)M'jj + sin2(/3 - a) M^] (sin^ aUjj + cos^ a Ml) , 

(32) 

where the last step was only possible because we are working under the assumption that the 
solution III is a minimum - therefore, M| > and Mfj > 0, and the inequality in eqs. (j32[l is 
not changed. Developing this inequality leads to a straightforward conclusion: 

Vi < Vm cos(2/3) < -1 . (33) 

A similar impossibility is found if one investigates the case Vji < Vm when III is a minimum. 
Which means that, if the solution III is a minimum, then it is certainly the global minimum 
of the theory. Then, the 2IIDM potential cannot tunnel from a phenomenologically acceptable 
minimum III, where all quarks are massive, to a deeper solution I or II, where either the up or 
down quarks could be massless. 

As a curiosity, we may also have a situation where the deepest minimum is either of the 
form I or II. In that case, there are two observations to make: (a) the solution III is necessarily 
not a minimum; and (b), there is the possibility that both solutions I and II are simultaneously 
minima. According to eqs. (|25p and ()27p . the specific values of the parameters will determine 
which of the two solutions corresponds to the deepest minimum. 



5.2 Normal vacua in the potential Vb 

We remind the readers that for the potential Vb we cannot choose a field basis so that 03 and 
{^13 , ^23} are simultaneously zero. The stationarity conditions for normal vevs in the potential 
Vb are very similar to those of eq. (I14p . namely 

2aivi + 03^2 + 4:biivf + 2(612 + ^33)'yi^^2 = 

2a2V2 + a^vi + Ab22V2 + 2{bi2 + b3s)vlv2 =0 . (34) 
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dV 
dV 



The presence of the 03 terms makes it impossible to solve analytically these equations. We can 
however follow a similar strategy to the one we used to treat the most general 2HDM potential: 
with polar coordinates v and /3, we obtain an equation for in terms of /3, 

2 _ 2ai +03 tan/? 

" cos2/3 2(612 + 633) tan2/3 + 46ii ' ^ ' 

and a quartic equation for tan /3, 

a3 622 tan'' /3 + [2 ai 622 - 02 (^12 + &33)] tan^ jS + [oi (612 + 633) - 2 02 fen] tan /? - 03 611 = . 

(36) 

This is an equation very similar to the one we studied in section [3l so the conclusions we 
reached there are still valid: this potential can have, at most, two pairs of non-degenerate minima 
away from the origin. 

To verify the previous statement, we performed a scan of the parameter space of the Vb 
potential searching for such minima. As before, our procedure was to randomly generate a 
partial set of parameters for the potential and two pairs of vevs - {vi , ^2}, such that v\ + V2 = 
(246GeV)2 for the A'^i minimum, and a second pair {vi , V2}, such that vf + has a value 
between (1 GeV)^ and (1000 GeV)^ for iV2. We then use the stationarity conditions ([M|) to 
determine, solving a set of linear equations, the parameters {ai , 02 , 03 , 633}. So that minima of 
the type A'^i have some phenomenological relevance we excluded the combinations of parameters 
which produced scalar masses too high (above 1000 GeV) or too low (below about 100 GeV). We 
found many such A'"i solutions, and for a small subset of those the minima A''i and A'^2 coexist. 
In fig. [3] we show the plot of the mass of the lightest neutral Higgs boson mass at the minimum 
A^i - the standard one, for which Mw = 80.5 GeV - versus the mass expected for the W boson 
at the normal minimum N2, when the deepest minimum found is N2. For many other points, 
we found that A'^i was the deepest minimum. The difference in depths of both minima can 
be large. The remarkable thing that this figure demonstrates is that we may be "living" in a 
perfectly reasonable A'^i minimum, where the known particle masses are those that have been 
measured; but "below" lies a deeper N2 minimum, with exactly the same gauge symmetries 
but very different particle mass spectra. In those circumstances, then, tunneling to the deeper 
minimum - and to a universe with very different particles - is, in principle, possible. We remark 
that there doesn't seem to be anything particular - any special combination of parameters of 
range of values of particle masses, for instance - characterizing either of the two regimes, deeper 
A'^i or deeper A''2. We have lightest Higgs masses ranging from 100 GeV to ~ 800 GeV for the 
A'^i minimum in either regime. 

There are two fulcral observations to make respecting these conclusions. Firstly, the per- 
centage of parameter space where both minima were found to coexist is extremely small: out 
of a generated 15 million A'^i minima, a A'^2 minimum existed in only about 37000, and of those 
in only about 1/3 of the cases was deeper than A'^i. Therefore, this feature of the 2IIDM 
potential is a rare occurrence, one such case for every thousand trials, but it exists nevertheless. 
Secondly, even if the A^2 minimum exists and is below A'^i, it is not clear whether that means 
the A^^i minimum is unacceptable or even "dangerous" . If the tunneling time from A^i to A'^2 is 
vastly superior to the current age of the universe, then A'^i would be an acceptable vacuum of 
the theory. The calculation of tunneling probabilities in models with more than one scalar is a 
very complex undertaking [T7], and outside the scope of the present work. 
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Figure 3: Lightest neutral Higgs scalar mass at the highest Ni minimum versus the mass of the 
W gauge boson at the deepest N2 minimum. 

6 Conclusions 

We have performed a thorough analysis of the structure of neutral minima in 2HDM potentials. 
We have shown that it is possible to have coexisting neutral minima in such potentials. Prom a 
careful study of the stationarity conditions we were able to establish the maximum number of 
possible minima that might exist in those conditions. For the 2HDM potential with explicit CP 
breaking, we concluded that there might only be one normal minimum with real vevs. However, 
for this potential, there is no physical distinction between stationary points with real vevs or 
with vevs which have a relative complex phase. We were then able to scan the parameter space 
of the potential and discover many combinations of parameters for which these two types of 
stationary points are indeed minima, with different depths. We chose a particulary curious case 
in which both of those minima predict identical values for the masses of the known particles 
(gauge bosons and fermions) but have completely different scalar spectra. 

In 2HDM potentials with spontaneous CP breaking we found that CP breaking minima and 
normal minima can coexist side by side in the potential, but the combinations of parameters 
corresponding to this situation are extremely rare. Further, requiring that the CP minimum 
describe the current known particle masses would imply a spectrum of scalars with very low 
masses, thus seemingly ruled out by experiments. 

For potentials where CP breaking is not possible, we can still have normal minima coexisting 
side by side. For the class of potentials Va, which have a Z2 symmetry, we found that it is possible 
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to demonstrate, analytically, that if there is a minimum where both fields $ have non-zero vevs, 
that minimum is unique and certainly the global minimum of the theory. However, for another 
class of potentials - Vb, which have a U{1) symmetry - we may have combinations of potential 
parameters for which two normal minima exist, with different vevs and different scalar spectra. 
In particular, we showed that it is possible to have the least deep of those minima with v = 246 
GeV, and the deeper one to have a completely different value for the squared sum of the vevs. 
This raises the possibility of tunneling between both minima. In fact, the results shown in 
this paper raise several interesting questions regarding cosmology: is it possible that tunneling 
between two normal minima, such as were found for the potential Vb, occurs in less than the age 
of the universe? Could we have a cosmological evolution described by a 2IIDM potential with 
spontaneous CP breaking in which the universe first "rests" in a CP breaking minimum, before 
tunneling or sliding down to the normal minimum that it currently occupies? And what would 
be the consequences of such an evolution in what regards questions of baryogenesis and matter- 
antimatter asymmetry? The study of tunneling in models with several scalars is complex, so 
these questions lie beyond the scope of this present work. 
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A Difference between the depth of the 2HDM potential at two 
normal minima 



As was mentioned in sec. [3l we have the possibility, in 2IIDM potentials, of having multiple 
stationary points with vevs of the "normal" type of eq. ([8j). Let us then consider two such 
stationary points, one with vevs {^i , ^2}, which we will call A^i, and another, N2, with vevs 
{vi,V2}- Let us also introduce the notation used in [4] , and define a vector X = [xi X2 X3 X4], 
containing the values of the x variables of eq. ([2]) at each of the stationary points. We also define 
the quantities V- = dV/dxi. For instance, Vi = ai -|- 2 5iiXi -|- 612X3 -|- 613X3 -|- 614X4. The 
vector y is thus defined as = [V-[ V2 V3 ^4], evaluated at each of the stationary points. The 
stationarity conditions of the 2HDM potential (eqs. (|14p to (jlSp for the most general potential, 
or the equivalent ones for the more restricted potentials) imply that, at the A^i stationary point 
we have 



vl 




V 



2vi V2 



Ni 



2vi V2 





(37) 



with analogous definitions at the A'^2 stationary point for the vectors X2 and V', with the obvious 
replacements Vi — > Vi. The quantity — V^/2vi V2, as was shown in ref. is related to the value 
of the squared charged scalar mass. Namely, we have 



2vi V2 



Ni 



vl + vl 



(38) 



with an analogous result for N2. Finally, let us use two more definitions introduced in ref. [1], 
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With this vector and matrix of parameters, we may write the value of the potential at each of 
the stationary points in a very concise manner. As was shown in [1] the value of the potential 
at A'^i is given by 



Vn, 



2 



- XI B Xi 
2 ^ 



(40) 



and the vector V at A'^i has a very simple expression: V = A + B Xi. Likewise, the value of 
the potential at N2 will be given by = \ X2 and V , at this stationary point, is given 
hyV' = A + BX2. 



With all necessary definitions introduced we may now demonstrate eq. ()19p . From the 
definitions of X and V and their values at the stationary points given by eq. (|37p . it is simple 
to see that we have 



X^V' 
XjV 



Xf{A + BX2) 
+ BXi) 



XlA + XfBX2 = 2VN,+XfBX2 
X^A + xJbXi = 2Vn2+xJbXi 



(41) 



where the last equality follows from eq. (j40p . and its analogue for the stationary point. Notice 
that, because the matrix B is symmetric, the two terms Xf B X2 and X2 B Xi are identical. 
So that, subtracting these two equations, we obtain 



— '^Ni 



l{x^v'-xfv') . 



Now, using equations ([37|) and ([38|) we can write 



(42) 



X^V 



[vl vl V1V2 0] ^ 



2viV2 





(vi V2 



V2 Vl] 



(43) 



and an equation entirely analogous to this one for Xj^ V' . From these two equations one obtains 
the result expressed in equation (fT9|) . 



B Classification of the several 2HDM potentials 



As explained in section [21 there are many different types of 2HDM potentials, depending on 
whether CP is or is not conserved, and on the types of symmetries that one imposes on the 
models. There is however a simple way of grouping those several potentials in specific categories, 
characterized by a single number, which we call the "index" of the potential. With the definitions 
of the real vectors X and A and the matrix B in the previous appendix, it is trivial to see that 
the most general 2HDM potential is written as 



V 



A^ X + - X'^ BX 
2 



(44) 



Now, B being a real and symmetric matrix, it can be diagonalised by a given orthogonal trans- 
formation O, such that 

"Si 
62 
53 
64 



OBO' 



(45) 
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Accordingly, the vectors X and A are transformed by the matrix O, 

X ^ X = OX , A A = OA , (46) 
so that the potential is now written as 

V = a-ix^ + hixj (47) 

with a sum on the index i assumed. Because i? is a 4 x 4 matrix and A and X vectors with four 
elements, we conclude that through the transformation O we can write the potential in terms 
of only eight quantities a and b - we say that this potential has index eight. Notice that O is not 
a basis transformation, and that we have not reduced the number of independent parameters of 
the potential - were we to study the stationary points of V, we would need the original 14 real 
parameters of the potential (or 11, with a suitable basis transformation) to do so. Equation (j47|) 
is nothing more than a simpler way of writing the potential. 

For the potential with explicit CP conservation, the matrix B has zeros on its fourth row 
and column, except the diagonal element, 644. Also, 04 = 0, which means that the last entry 
of the vector A is zero. In this case, the transformation O that diagonalises B has zeros in its 
fourth row and column, except the (4 , 4) element, which is equal to 1. Then, the rotated vector 
A = O A still has a zero in its fourth entry. The rotated matrix B has a priori four independent 
eigenvalues, so this potential ends up being written in terms of seven parameters - three and 
four 6j. The potential with explicit CP conservation has therefore index seven. 

What about the potentials Va, Vb-, Vc-, for which extra symmetries have been imposed? 
Well, for Va and Vb the matrix B has further zeros, since 613 = 623 = 0. The matrix O is 
therefore block diagonal, with a 2 x 2 matrix in its first two rows and columns and the identity in 
the third and fourth positions. The transformation O therefore does not affect the values of 03, 
^4; ^33 and 644. The transformed elements therefore satisfy the conditions 03 = 03, 04 = (04 
is already zero for these potentials, since they explicitly preserve CP), 63 = 2 633 and 64 = 2 644. 

The potential Va (eq. ([5])) has 03 = and 633 / 644. Therefore, after the transformation 
O, it will have index six - two Oj and four hi parameters. If one includes a soft breaking term 
(an 03 term) then there is an extra parameter - the softly broken Va has index seven. For the 
potential Vb, because 633 = 644, two of the eigenvalues of B are equal. Therefore, there are 
three different Oj parameters and three bi ones - Vb has index six. Notice that, had we allowed 
a massless axion in this model (or in Va) , the index number would have been reduced by one. 
Finally, for the potential Vc of eq. ([7]), the diagonalisation of the matrix B and corresponding 
rotation of the vector A leads to the conclusion that this model has index six, and that the 
softly broken Vc has index seven. The study of the stationary points of the potential through 
the diagonalization of an analogue of matrix B was done in ref. ^ . An interesting observation 
is that all the potentials where CP can be broken (explicitly or spontaneously) have index larger 
or equal to seven. 
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